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ABSTRACT 

We describe a new iteration method to estimate asteroid coordinates, which is based on the 
subpixel Gaussian model of a discrete object image. The method operates by continuous pa¬ 
rameters (asteroid coordinates) in a discrete observational space (the set of pixels potential) 
of the CCD frame. In this model, a kind of the coordinate distribution of the photons hitting 
a pixel of the CCD frame is known a priori, while the associated parameters are determined 
from a real digital object image. The developed method, being more flexible in adapting to 
any form of the object image, has a high measurement accuracy along with a low calculating 
complexity due to a maximum likelihood procedure, which is implemented to obtain the best 
fit instead of a least-squares method and Levenberg-Marquardt algorithm for the minimisation 
of the quadratic form. 

Since 2010, the method was tested as the basis of our CoLiTec (Collection Light Tech¬ 
nology) software, which has been installed at several observatories of the world with the 
aim of automatic discoveries of asteroids and comets on a set of CCD frames. As the result, 
four comets (C/2010 XI (Elenin), P/2011 NOl(Elenin), C/2012 SI (ISON), and P/2013 V3 
(Nevski)) as well as more than 1500 small Solar System bodies (including five NEOs, 21 Tro¬ 
jan asteroids of Jupiter, and one Centaur object) were discovered. We discuss these results that 
allowed us to compare the accuracy parameters of a new method and confirm its efficiency. 

In 2014, the CoLiTec software was recommended to all members of the Gaia-EUN-SSO 
network for analysing observations as a tool to detect faint moving objects in frames. 

Key words: methods: data analysis; minor planets, asteroids, comets; techniques: image 
processing. 


1 INTRODUCTION 

There are many methods for determining an asteroid position dur¬ 
ing observations with a CCD camera. For example, the Ful l Width 
at Half Magnitude (FWHM) approach iOarv & Healvl2003i . which 
is based on the analytical description of the object images on the 
CCD frame, as well as other methods, in which the position of an 
object’s maximum brightne ss on a CCD image is taken as its coor¬ 
dinates jMiura et alj EoOSh l. 

Most of these methods have a common feature. They use 
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PSF-fitting (Point-Spread Function) to approximate the object im¬ 
age and get the information about regularities in the distribu- 
tion of the regis t ered p h otons on CCD fr a me(see, for detail. 


tion ot me regis t ered p h otons on V .r .1-) tr a me(see, tor detail, 
lYanagisa wa et alJ bOOSl). iG ural et a ll (l2Q05h. IVeres and Jedicke 


201^, Dell'Pro et alj I ^12h . Lafreniere et^l EoO'E ' IZachariaa 

1 2010h ). Among the models of photons distribution, which are 


used more often, we not e the two-dimensio n al Gaussian model 
I Veres and Jedickel (|201^. Zachariaj 2010h. Jogesh Babu et all 


2008h ). Moffat model teauej 1 20K h 


Izmailov et al. 


koioh ) 


Lorentz model jZachariasI i2010t) . lizmailov et alj ^20101) ). These 
models are usually described by the continuous functions, while the 
CCD images are discrete ones. Such an approach was criticised rea- 


© 2014 RAS 
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sonably bv lBauej l l2009l) . The principal disadvantage is that these 
models work well only with a large amount of the data. It leads to 
the fact that, firstly, the computation process becomes much more 
complicated, and, secondly, the problem related with the adequacy 
of the used estimations of PSF parameters cannot he solved. As 
a result, the error of the coordinate determination of the observed 
celestial objects has to he increasing. 

In addition to the aforementioned disadvantage, the existing 
methods do not pay sufficient attention for taking the noise com¬ 
ponent of the object image into account. It is assumed that its reg¬ 
istration and compensation are performed dur ing the preliminary 
stage of image processing dOural et al.l ( l2005h l or that the object 
image is exem pted from noise accord i ng to the accepted si gnal- 
to-noise model jLafreniere et alJbOOTll . lizmailov et'^ ( l2010h l. At 
the same time, the error introduced by the operation of removing 
the noise component from the object image is not considered in the 
subsequent procedure of coordinate determination. 

The errors of CCD astrometry are traditionally divided into the 
instrument, reduction, reference catalogue and measurement errors. 

The first ones traditionally include errors of instrumental pa¬ 
rameters as shutter delay and clock correction, which result in in¬ 
correct timing. The second error category (reduction) is associated 
with the method related to the standard and measured coordinates 
and depends on the choice of the algorithm relation between these 
coordinates. For example, being not sufficient for the wide-angle 
astrographs such an error depends strongly on the choice of the type 
and degree of the polynomial approximation in Turner method. It 
is important to note that a systematic error of timing is not shared 
with the coordinate error along the tracking an object in the sky, so 
it must be caught clearly with a reliable shutter sync. 

The reference catalogue errors are divided into three main 
classes: zonal errors (systematic errors of the reference catalogue); 
coordinate errors of reference stars in the catalogue epoch; proper 
motion errors of reference stars. Therefore, a choice of the refer¬ 
ence catalogue is very important. For example, Hipparcos or Tycho 
catalogues had no errors in the epoch of 1990.0, because an in¬ 
trinsic accuracy was at the level of millisecond of arc, but for the 
present epoch there is a necessity to take into account the proper 
motion errors of reference stars. The solution of this problem, i.e., 
creation of new huge database of proper motions of stars, is one of 
the task of the GAIA mission. It is worth noting that the choice of 
reference catalogue is not so important, when monitoring observa¬ 
tions of sky are conducted with the aim of discovery of new Solar 
system small bodies, because an intrinsic accuracy of a catalogue 
should not be necessarily maximised as compared with those for 
the following tracking the discovered object. 

The measurement error is related, first of all, to the determina¬ 
tion of coordinates of the image centre or the accuracy of a digital 
approximation of the CCD image (fitting). The attempt to improve 
the fitting may not lead to expected results if the reference cata¬ 
logue errors and timing are not taken attentively into account. 

Each of the aforementioned factors could be the source of both 
systematic and random errors. Any attempt to reduce one of these 
errors is impossible without control of other factors. Therefore, the 
task of the observers is to be responsible for monitoring all the pos¬ 
sible sources of errors. The aim of this paper is to help the observers 
to refine the coordinate measurements of the object image on the 
CCD frame and to control the errors of the measured coordinates. 

With this aim, we developed a new method for accurately esti¬ 
mating asteroid coordinates on a set of CCD frames, which is based 
on the subpixel Gaussian model of the discrete image of an object. 
In this model, a kind of the coordinate distribution of the photons 


hitting a pixel on the CCD frame is known a priori, while its param¬ 
eters can be determined from the real digital image of the object. 
Our method has low computational complexity due to the use of 
equations of maximum likelihood, as well as the proposed model is 
more flexible, adapting to any form of the real image of the object. 
For example, in fact PSF is a super-pixel function because it de¬ 
scribes the changes in the brightness of a pixels of celestial object 
image. We propose to use the density function of coordinates of 
hitting photons from celestial object instead of the PSF. This func¬ 
tion is sub-pixel one. To get the Point-Spread Function from this 
function, it should be integrated over the area of determination of 
each image pixel of object or compact group of objects. It turns out 
that sub-pixel models are more flexible and may describe the real 
image more adequately. This effect does not occur in case of bright 
objects. But applying a more flexible model for the faint objects, we 
are able to improve the accuracy of the measurements (for exam¬ 
ple, as comparing with Astrometrica) by 30-50% (more details are 
given in the discussion). Our method also takes into account the 
principal peculiarities of the object image formation on the CCD 
frame along with the possible irregular distribution of the residual 
noise component both on the object image and in its vicinity. 


Some generalisatio ns of the proposed method are presente d 


in our previous works JSavanevvchl 1 1999 

). Savanevvchl (2006b. 

ISavanevvch et alj (|201(]I). Savanevvch et al 

(I 2 OI 1) Vavilova et al.l 

(2011). S avanevvch et al. 

i2012h. Vavilova et all (l2012ta)h. 

IVavilova et al. (2012(b)h. Savanevvch etalj i2014|)). Since 2010. 


due to its application with the CoLiTec (Collection Light Tech¬ 
nology) software, which was installed at several observatories of 
the world, four comets (C/2010 XI (Elenin), P/2011 NOl(Elenin), 
C/2012 SI (ISON), and P/2013 V3 (Nevski)) and more than 
1500 small Solar System bodies (including five NEOs, 21 Trojan 
asteroids of Jupiter, and one Centaur object) were discovered. 
These results confirm the efficiency of the proposed method. The 
main stages of image processing with the CoLiTec software are 
presented in Fig.[T] 

The structure of this article is as follows: we describe a prob¬ 
lem statement and the method in Chapters 2-3, respectively. The re¬ 
sults based on the testing of this method with the CoLiTec software 
are presented in Chapter 4. We compare and discuss the accuracy 
and other parameters for determining the position of the faint ce¬ 
lestial object on the CCD frame obtained by the proposed method 
and others in Chapter 5. Conclusive remarks are given in Chapter 
6 . 


2 PROBLEM STATEMENT 

If the exposure time is small, the shift of asteroid position in the 
sky can be ignored. In this case, the asteroid and field stars are 
imaged as blur spots rather than points on the CCD frame. It is 
postulated that the coordinates of the signal photons from asteroids 
and stars hitting the CCD frame have a circular normal distribu¬ 
tion with mathematical expectations xt, yt and mean-square error 
(MSE) (jph- 

It is supposed that a preliminary detection of asteroid has al¬ 
ready been conducted before the determination of its coordinates. 
The result of this detection is a preliminary estimation of the aster¬ 
oid position on the CCD frame, namely, the determination of coor¬ 
dinates of the pixel, which corresponds to the maximum brightness 
peak on the asteroid image. We name the set of pixels around this 
pixel as the area of intraframe processing (AIFP). Thus, the AIFP 
size {NiPS, in pixels) is much larger than the image of the asteroid. 
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New method for estimation of asteroid coordinates 3 



Figure 1. Main steps of object image processing in the CoLiTec software: a) exclusion of stationary objects; b) detection of moving objects; c) analysis of 
moving objects, where Al, A2, and A3 outline coordinate deviations of moving objects from their trajectory. 


The original CCD image of the celestial object contains harm¬ 
ful interferences such as the read noise, dark currents, irregu- 
larit y in the pixel-to-pixel sens i tivity, sky bac kground radiation, 
etc. l lFaraii and MacLeanI i2006h . lHarrij l ll99(]h l. Hence, the CCD 
frame, can be represented as an additive mix of the images of celes¬ 
tial objects and a component, which is formed by this generalised 
interference. Within the scope of the whole CCD frame the inter¬ 
ference component has a complex structure. However, in a small 
vicinity of the studied asteroid image, such interference component 
can be described with a good accuracy as a plane with an arbitrary 
slope. Such a representation describes well the interference com¬ 
ponent especially if there is a bright object near the vicinity of the 
studied AIFP. 

The output signals of the CCD matrix pixels Nips are easily 
reduced to the relative frequencies of the photons hitting the 
pixel on the frame: 


^ikt 


-^ikt 


Nips 

Aikt 


i,k 


( 1 ) 


where Aikt is the brightness of the ik pixel of the CCD matrix. 
Then, the result of the observation is the set U = 
..., i^Nipst of tho relative frequencies, which are in¬ 
dependent of each other. The theoretical analogues of the measured 
relative frequencies are the probabilities iiikt{Q) that during the ex¬ 
posure time the photons hit the ifc*^-pixel of the CCD matrix with 
the borders Xbegi, Xendi in the coordinate x and ytegk, Vendk in the 
coordinate y on the frame. It is supposed that the angular sizes 
of the pixel, Aa; and Ay, are the same in both coordinates x and y. 

Thus, the problem statement is as follows: it is needed to de¬ 
velop a method of maximum likelihood estimation of asteroid co¬ 
ordinates on the CCD frame using the set of relative frequencies 
v*kf It is believed that the likelihood function is differentiable in 
the vicinity of its global maximum, and its initial approximation is 
also in the same vicinity. The set of the estimated parameters 0 in¬ 
cludes the asteroid coordinates xt and yt on the d^ frame and the 
mean square error of the coordinates of photons hitting the CCD 
matrix, aph- 


To introduce a new method, we use two functions. The den¬ 
sity distribution of a normally distributed random variable a with 
mathematical expectation rriz and dispersion is determined by 
the expression: 


Ndmz,ix‘^) =-^=exp{--^{z ( 2 ) 

The probability that a random variable z is within the 
closed interval [zb<,g , Zi,nd] is: 


Fzi(m^, = 


f 

J 


N2{mz-,a)dz. 


(3) 


3 TASK SOLUTION 

For achieving the maximum accuracy of the estimations of the ob¬ 
ject position on the frame, the discretisation factor needs to be taken 
into account, because we should estimate the continuous parame¬ 
ters (coordinates of objects) at the discrete set of the measured val¬ 
ues (the brightness of the CCD matrix pixels). The general view of 
the maximum likelihood estimation of the object position can be 
expressed by: 


^ikt dviktjQ) ^ , 

ViktiQ) dQm ' ’ 

where 0 is the set of the estimated parameters Xt,yt, and aph- 
The relation between the probability i/ikt (0) that photons hit 
the pixel (Eq. 4) and the function of coordinates distribution 
f{x,y) of the incidence of photons from the object on the CCD 
matrix has the form: 




r^endi P 


Vendk 

Vbegk 


f{x,y)dxdy. 


(5) 


After the compensation of the noise component on the CCD 
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image, the function f{x, y) could be expressed as the weighted mix 
of normal and uniform probability distributions: 

/(a:,t/,e)=po + ^ 


^'^ph 

where pi = 1 — po is the relative weight of the signal photons 
of the object; po (0 ^ Po < 1) is the relative weight of the residual 
noise photons of the CCD matrix after the compensation of the 
flat generalised interference; xt and yt are the object coordinates 
on the frame at the time corresponding to the mathematical 
expectations of the coordinates of incidence of the signal photons. 

The probability (Eq.5) that photons hit the CCD matrix pixels 
can be written as follows: 


— ^iktnoise likts^ (7) 

where likta = PiFa:i{xi\alf^i)Fyk{yt, <xlh) is the probability 
that signal photons Fkt hit the pixel of the CCD matrix; and 
liktnoise = A^copO i® the probability that the noise residuary 
photons hit the pixel of the CCD matrix; Accd ~ = Ay. 

The derivative from the probability (Eq. 7) in the x coordinate 
is determined by the expression: 


dvikt{Q) 

dxt 

where 


PiFyk{^yu'Xph)Fj;i{xt\(Jy^') 

- 


^ph 


-Xt), (8) 


tU-xz rrix T _ . Q xendt {jXlx , tT j Axbegi {xklx , O j j 

k xi ybrix ; rr j 

is the local (on the closed interval [xbegr, xendi]) mathemat¬ 
ical expectation of the normally distributed random value x with a 
mathematic expectation rrix and dispersion cr^. The derivative from 
the probability (Eq. 7) in the y coordinate has the same expression 
as Eq. 8. 

The system of equations of the maximum likelihood for the 
studied AIEP pixels in the case while the asteroid position is esti¬ 
mated only, will take the form: 


^ IKt XI ^ 

“ ^^IPSs * X ■ 

^ i,k ^ikt 
y^^IPSs * \ loc 

^i,k ^ikt^ikt'^yk 


(9) 


yt 


^i,k 


'ikt 


^ikt 


where AjpSg is the quantity of pixels in the part of AIEP area 
(the area, where the signal from the object is expected); xt and yt 
are the estimations of asteroid coordinates; and Xtkt is the part of 
photons from the celestial object in the ik*^ pixel of the CCD 
frame. The latter value is determined by: 


Xikt — 


PlFykitjt', <Xph)Fxi{{Xt’, tXph) 
t^ikti.^) 


( 10 ) 


To estimate the MSE of the coordinates of the photons hit¬ 
ting the CCD frame from the asteroid, we use the equation of the 
maximal likelihood: 


- 2 
^ph 


E ^IPSs 

i,k 


v*ikAikt{{m}xi - xt)^ + {rriyk 


2e: 


NiPSb 


^ikt^ikt 


yt)^) 


( 11 ) 


We cannot exclude completely the generalised noise interfer¬ 
ence. By this reason, to take into account the relative weight of the 
signal photons, we use a standard estimatio n of the weights of the 
discrete mix of the probability distributions (ILo etalJdTOOlh ): 


Pi 


1 

AlpSa 


^IPS s 

Mkt ; po = 1 — pi ■ 
i,k 


( 12 ) 


Therefore, the local mathematical expectation of coordinates 
of the object position is a function of the relevant coordinates, 
and Eq. 9 gives a system of transcendental equations that can 
be solved by the method of successive approximations (see, e.g., 
iBurden and EairesI d201(i) ). 

The algorithm of the estimation of object coordinates consists 
of two successive operations. 

The first operation is to split the statistics of the AIEP pixels 
into the statistics of the signal and statistics of the residual inter¬ 
ference. It is performed for pixels from the AIEP area where the 
object is expected. According to the 0 values calculated from the 
previous iteration, the photons of pixel are divided into those be¬ 
longing to the object and to the residual interference. The photons 
belonging to the object are analysed to determine estimation of its 
position. Thus, coordinates of the local maximum in the object im¬ 
age, around which the AIEP area is formed, are used as the initial 
approximation. The result of this operation is a set of split coeffi¬ 
cients Aifcf 

The second operation provides estimation of the object coor¬ 
dinates based on the statistics, which are obtained during the opera¬ 
tion of splitting. It is conducted in a strongly determined way from 
Eq. 9 to Eq. 12. The result 0„ of this operation serves as an ini¬ 
tial approximation for the operation of splitting at the next iteration 
step. The iteration process is continued until the difference between 
0n and 0n_i becomes smaller than the predetermined value, for 
example, 0.1% of the angle size of a pixel. 

The analysis of the iteration process shows that its conver¬ 
gence is provided while the following conditions are fulfilled: 


dx = <6,dy= < 6, (13) 

(7x O' y 

where dx and dy are the relative distances between the initial 
and actual positions of the object; xo and t/o are the initial approx¬ 
imation of the object coordinates, Xtme. and ytma are the actual 
object coordinates; and ax and Oy are the MSEs of coordinates of 
the signal photons hitting the CCD matrix. 

The observations based on the proposed method have shown 
that this condition is almost always fulfilled for the real images of 
asteroids and stars, and, in most cases, the relative distance is not 
more than 1.0 -T 15. 

The opportunity to divide the AIEP area into the interference 
area (pixels that have registered photons only from interferences) 
and the object area (pixels that have registered photons from the ob¬ 
ject and interference) yields a more simple and reliable algorithm 
of estimation of the flat interference component (as compared with 
the estimation in the common system of maximum likelihood equa¬ 
tions). Namely, there is an independent estimation by the method 
of least squares (MLS). Thus, the density of the coordinate distri- 
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bution of the photons from the residual interference will represent 
an equation of the plane with an arbitrary slope: 


fnoisei^^^y') — ^noise^ Bnoisey “t“ Cn 


(14) 


The probahility that these photons will hit the pixel can 
be given by analogy with Eq. 5: 


^iktn 


,(e„ 


\ - 4 ^^^. rp.j _ 1 _ fj.-j I 

e) — -^no^se'^^K i ^no^sey^K \ '^nois 


(15) 


where I'ii^tnoise is th® measured frequency of the noise 
photons hitting the pixel of the CCD matrix; — 

NfnnAnoise^ — NnnnBn 


-yint 

''noise 




= {A^ 


^CCD ^noise > 

are the integral parameters 


of the flat noise component and its vectors; Xikt ~ 

yikt = are the average values of coordinates of the 

ik^'^ pixel. 

Thus, the probabilities that noise photons will hit the pixels 
of the studied AIFP depend linearly on the angle coordinates of 
the centres Xj and yj of these pixels and represent, by themselves, 
a plane with the integral parameters and C^Z\se- 

It is worth noting that the pixels containing the supposed object 
image should he eliminated before the determination of the noise 
parameters. 

The integral parameters of the fiat interference component 
A'iifise, and Clifisfi can be determined with a linear MLS 

estimation: 


where 




^noiae 

= (F^F)- 

■ 

J- L^nozse? 

(16) 

Xl 

xi 

^Nxpsnoise 


yi 

- yi 

y^IPSnoise 

(17) 

1 

... 1 

1 



where Xj and yj are the angular coordinates of the pixel, 
which is used to estimate parameters of the fiat interference com¬ 
ponent; Ih® number of AIFP pixels, which do not 

contain the object image. 

To obtain the integral parameters (Eq. 16 - Eq. 17), only the 
AIFP pixels not belonging to the areas with the object images 
(Nipsnoiae ^ {Njps — Nipsa)) should be used. To exclude the 
influence of the anomalous emissions of the brightness in the pix¬ 
els, we apply two iterations of MLS. In the second iteration, we use 
only those pixels, for which the value of the obtained frequency 
fiktnoiae Satisfies the following condition: 


l^iktnoiae ^iktno 


€ 


^ kn 




2-^i,k 


'-ft, 


Nips„ 


(18) 


where k„oiae is the threshold coefficient for removing the pix¬ 
els, which do not satisfy this condition, for example, knoiae='i', 

2 

- are the standard deviations 
of the flat interference component; and Nktrioiae = + 




°)(L*fctnoise- 


Nipsnoiae 


Bttoiaeykt IS the Smoothed estimation of the measured fre¬ 

quency of the ikt'^ pixel that is a part of the Nipsnoiae value given 
in Eq. 16. 

The number of pixels Nipsnoiae determining the set Unoiae 
is reduced on its quantity, for which the condition (Eq. 18) is not 
satisfying. The process is repeated until one of the following con¬ 
ditions is completed: 1) the module of difference of the two related 
values of standard deviations becomes less than a certain value; 2) 
the number of pixels Njpsnoiaa becomes less than a given number; 
3) the number of iterations exceeds a predetermined limit. 

The next step is as follows: the obtained values of the integral 
parameters of the flat interference component are subtracted from 
the object signal of the given AIFP: 


* _ * / ^int ^int . ^int \ /-t q \ 

^ikts ^ikts noise Xits + ^noiseVkts ' ^noise)^ 

where i^ikts is the measured frequency of the photons hitting 
the ik*^ pixel from the AIFP area; xua and ykta are the angular 
coordinates of the ik^^ pixel from the AIFP area. 

The AIFP area for the calculation of the flat interference com¬ 
ponent of the object image is set by the operator (the default 31x31 
pixels). If the size of the object image is larger for this area, then 
the size of the area for the flat interference component is taken as 
the size of the object image multiplied by two (it is also set by the 
operator). As for the area of frame for fitting, we note that the fit¬ 
ting is carried out on the pixels that belong to the object image. 
Determination of these pixels is conducted through the delineation 
procedure, i.e. the fitting is carried out not for the fixed area but for 
the area which depends on the size of the object image. 

The principal stages of the object image processing in the Co- 
LiTec software are demonstrated at the pipeline in Fig. 1. They 
include: a) exclusion of stationary objects; b) detection of mov¬ 
ing objects; c) analysis of moving objects, where Al, A2 and A3 
outline the coordinate deviations of moving objects from their tra¬ 
jectory. If it is necessary, the measurements can be stacked, but 
only un-stacked frames are statistically processed. Image service 
files (flat - bias - dark) can be used during the calibration process, 
but we offer an alignment frame option, which is commonly used. 
For example, to align frames, we used a high-pass Fourier filter in 
the earlier versions of the CoLiTec software. Currently we are us¬ 
ing a median filtering, which has practically the same quality but 
it is substantially faster. To select the moving objects, the measure¬ 
ments (blobs) are formed in all the selected object images. After 
this, frames should be identified with each other and coordinates of 
all the measurements should be led to the basic frame. 

The algorithm of the proposed method, which could be helpful 
to implement it, works as follows: 

1. To form a studied square area of infraframe processing 
(AIFP) with a side of I pixels (NiPSs=s^ ) and a square region 
of the presupposed existence of images of celestial objects with a 
side of s pixels (s << I ). The centres of these regions are the local 
maxima of the image of the object (asteroid, comet, etc.) discovered 
previously. 

2. To conduct a multipass MLS-estimation of parameters 
of the interface noise component according to Eqs. 15 and 17. 
Wherein, only those pixels are processed at the next MLS- 
estimation pass, for which the value of the relative frequency 
(Eq. 1) satisfies the condition of Eq. 18. This do-while loop is re¬ 
peated as long as the module of difference of the two standard de¬ 
viation values obtained sequentially becomes smaller than a prede¬ 
termined value. 

3. To exclude the noise photons from the potentials of pixel 
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in a region of the presupposed existence of images of celestial ob¬ 
jects according to Eq. 19 (MLS-estimations of parameters of the 
interface noise component). 

4. To estimate the coordinates of the object position according 
to Eq. 19 on the digital image, from which the noise photons were 
excluded earlier. Initial values of the coordinates of object position 
should be equal to the coordinates of the central local maximum of 
the AIFP image (AIFP centre): 

4.1. to calculate the coefficients of splitting \ikt according to 
Eq. 10; 

4.2. to estimate the coordinates of object position on the CCD 
frame according to Eq. 9; 

4.3. to get the standard deviation estimate of the photons hit¬ 
ting from the object according to Eq. 11; 

4.4. to determine values of the weights of the discrete mix of 
the probability distributions according to Eq. 12; 

4.5. to compare the actual value and value obtained in the pre¬ 
vious step. If the difference between them is greater than a prede¬ 
termined value, the actual value goes to the item 4.1 as estimates 
of the previous step; otherwise, it goes to the algorithm output as a 
result of its work. 

The important step of the automated pipeline is also the cre¬ 
ation of a catalogue of stationary objects (blobs) in the frame series. 
All the objects of this catalogue should be excluded in next stages 
of the processing but a few objects should be left as reference stars 
for astrometry. Thus, the measurements, which are absent in the 
catalogue of stationary objects, have been used to select the moving 
objects. With this aim, the original method of collecting light (col¬ 
lection light technology) is provided. Candidates to the small Solar 
system bodies (e.g. asteroids) selected by this method should have 
about the same brightness on different frames, and their coordi¬ 
nates should not deviate significantly from their average trajectory. 
In case of a partial occultation of a star by an asteroid, more often 
they stand out as two different celestial objects, wherein the posi¬ 
tion accuracies of the star and asteroid are falling. In general, the 
detection of an asteroid is provided at the step of blobs, while the 
images themselves are not involved into processing (in this course, 
the blobs are formed according to the object images). 


4 RESULTS AND DISCUSSION 

The proposed method provides accurate estimation of asteroid co¬ 
ordinates on a set of the CCD frames and is the basis of the CoLiTec 
software, which was installed at several observatories of the world 
since March 1, 2009 (see, also, Fig.l). 

In April, 2010, the CoLiTec was installed at the Andrushivka 
Astronomical Observatory (A50), Ukraine, and already in May, 
2010, two asteroids were discovered (it was the first discovery of 
asteroids in the automatic mode at the observatories of the CIS 
countries). On November 27, 2010, this software was installed at 
the ISON-NM (H15, New Mexico, USA), and on December 10, 
2010, the comet C/2010 XI was discovered by L. Elenin. In July, 
2012, the CoLiTec software was installed at the ISON-Kislovodsk 
(DOO, Russian Federation), and on September 21, 2012, the comet 
C/2012 Sl(ISON) was discovered by V. Nevsky and A. Novi- 
chonok, the amateurs of astronomy. The digital images of discov¬ 
eries of these comets are given in Fig.|^ and Fig.|2j5, respectively. 
As we mentioned in the Introduction, in total 4 comets and more 
than 1500 asteroids have been discovered with this software since 
2010 (two examples of the asteroid discoveries are demonstrated in 
Fig. |2j; and Fig. [2ll, where the latter object discovered at the An¬ 


drushivka Astronomical Observatory on April 24, 2011, is marked 
as unusual K11H52Y asteroid in the MFC circular). 

Since this method was already being operated during the as¬ 
teroid observational surveys, we are able to provide a comparative 
analysis of the statistic parameters of accuracy estimations by this 
method and other results of thirty observatories that are most pro¬ 
ductive in terms of the number of asteroid observations in 2011- 
2013 (Tables 1-3). The observatories, which work with only one 
object in the centre of the CCD frame during the calm and suit¬ 
able phase of the Moon, are excluded from our analysis. The ob¬ 
servatories, where the CoLiTec software based on this method was 
installed, work during gusts and when it is calm. Moreover, the 
method is adaptive in such a way that there is no problem in auto¬ 
matic processing the object images on the CCD frames in the centre 
and on its edges as well as to provide measurements of many ob¬ 
jects on the single frame. 

Therefore, in the years 2011-2013, such obse r vatorie s 
as the ISON-NM Observatory (HI5) jElenin et al.l ( l2012h . 
lElenin et al.l (l20 13h). And r ushivk a Astronomical Ob se rvator y 
(A50) j lvashchenko et al.l (1201 ih . Ilvashchenko et al.l ||2012|) . 
Ilvashchenko et al. (2013)), and ISON-Kislovodsk Observatory 
(DOO) acted as the users of the CoLiTec software dElenin et al.l 
( I 2 OI 4 I) ). In the ranking of the most productive observatories 
worldwide in 2012, based on the number of discoveries of the 
small Solar System bodies, the users of the CoLiTec software 
had the 3rd, 13th, and 22nd places, respectively. As for the final 
report of 2011-2012, the ISON-NM Observatory (H15) holds the 
7th place both by the number of measurements and priority of 
discoveries. 

In Table 1 to 3, the total numbers of measurements and ob¬ 
jects (Column 3) as well as of the asteroid discoveries (Column 4) 
are given according to the circulars of 20 11, 20 1 2, and 2013 of the 
Minor Planet Center (MFC), respectively ( lMPd(l2013h ). The statis¬ 
tic parameters of these measurements are taken from the MFC site. 
The following parameters are also indicated in the tables for each 
observatory: diameter D of the primary mirror of the telescope, in 
meters (Column 5); scale Spix of the pixel image, in arc seconds 
(Column 6); the average residuals and Aj of object positions 
in the right ascension a and declination 5 at a predetermined time 
(Column 7); standard deviation estimations (Ja and as of object po¬ 
sitions in the right ascension a and declination <5 at a predetermined 
time (Column 8); standard deviation estimations a" of object po¬ 
sition (Eq. 20), in arc seconds (Column 9); standard deviation es¬ 
timations cTpia; of object position (Eq. 21), in pixels (Column 10); 
module of the Average Residuals of object position Measurements 
(Eq. 22), ARM (Column 11). To calculate some aforementioned 
parameters, we applied the formula as follows: 


a” = 0, 5(CTa -f as); 


( 20 ) 


^pix — 



( 21 ) 


ARM = y'jAcj^-bjAajL (22) 

The data analysis related to the accuracy parameters of object 
position for the most productive observatories in the number of as¬ 
teroids measurements in 2011-2013 is illustrated in Figure 3. 

The observatory-partners of the CoLiTec program hold lead¬ 
ing positions in their class of telescopes as related to the parameter 
of module of the average residuals of measurements (Figures 3a). 
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Figure 2. The most known small Solar System bodies, which were detected and discovered with the CoLiTec software based on the subpixel Gaussian model: 
a) the long-period Comet C/2010 XI (Elenin), having an image size of 5 pixels, has been shifted by 7 pixels on a set of four CCD frames; b) the Comet C/2012 
Si (ISON), having an image size of 5 pixels, has been shifted by 3 pixels on a set of four CCD frames; c) Centaur object 2013 ULIO, having an image size of 
5-6 pixels, has been shifted by 2 pixels; d) the K11H52Y asteroid, having an image size of 5 pixels, has been shifted by 45 pixels on a set of three CCD frames 
(it is described as unusual in the MFC circular www.minorplanetcenter.net/mpec/Kl 1/Kl lJ02.html) 


In 2011 and 2012, this parameter for H15 telenin et alj ( l2012[) ) and 
A50 observatories was equal to 0.06”. At the same time, these ob¬ 
servatories were not in the list of the best ones as related to the 
parameters of standard deviation estimations of object position (in 
arc seconds) (see. Figures 3c and 3d). In 2011, the values of stan¬ 
dard deviation estimations a" of object position for the mentioned 
observatories were equal to 0.515"(H15) and 0.51” (A50). In 2012, 
these values were equal to 0.515" (H15), 0.49” (DOO), and 0.48” 


(A50). The reason for such a deterioration of the results, in addi¬ 
tion to the size of the aperture of a telescope, is the pixel scale of the 
used CCD matrix. To take this factor into account, the observatory- 
partners of the CoLiTec program decided to consider the parameter 
of standard deviation estimations of object position, in pixels, 
for accurate estimation of asteroid coordinates on the CCD frame 
as a principal one during observations. 

The parameter of standard deviation estimations cTpia, of the 
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Table 1. The accuracy parameter of thirty observatories that were the most productive in the number of asteroids measurements in 2011 according to the MPC 
data. 


(1) 

N 

(2) 

Observatory code 

(3) 

Measurements, objects 

(4) 

Discoveries 

(5) 

D, m 

(6) 

Spix 

(7) 

R.A.(Aq,/(Tc,) 

(8) 

Decl./A^/o-j) 

(9) 

a" 

(10) 

^pix 

(11) 

ARM 

1 

G96 

2106367, 382737 

21770 

1.50 

1.00 

-0.01 +/- 0.32 

-0.04 +/- 0.28 

0.300 

0.30 

0.041 

2 

704 

1956368, 279129 

495 

1.00 

2.20 

0.25 +/- 0.66 

0.43 +/- 0.64 

0.650 

0.29 

0.497 

3 

F51 

1557902, 351923 

13628 

1.80 

0.30 

0.05-t/-0.16 

0.06+/-0.17 

0.165 

0.55 

0.078 

4 

703 

1512387, 259412 

2995 

0.68 

2.60 

-0.21 +/- 0.67 

0.17+/-0.68 

0.675 

0.25 

0.270 

5 

691 

811571, 154495 

8356 

0.90 

1.00 

-0.16+/-0.33 

0.10+/-0.30 

0.315 

0.315 

0.189 

6 

E12 

219903, 52808 

327 

0.50 

1.801 

-0.04 +/- 0.49 

0.32 +/- 0.48 

0.485 

0.26 

0.322 

7 

645 

208656, 45961 

7 

2.50 

0.39 






8 

D29 

185303,43414 

318 

1.04 

1.70 






9 

C51 

162900, 15412 

23 

0.40 

2.75 

0.06 +/- 0.57 

-0.03 +/- 0.65 

0.610 

0.22 

0.067 

10 

HIS 

154970,37495 

768 

0.45 

2.00 

-0.03 +/- 0.49 

0.06 +/- 0.54 

0.515 

0.25 

0.067 

11 

106 

75340, 18093 

73 

0.60 

2.00 

0.04 +/- 0.36 

-0.11 +/-0.35 

0.355 

0.17 

0.117 

12 

291 

70355, 19028 

646 

1.80 

0.60 

-0.13+/-0.36 

0.15+/-0.27 

0.315 

0.52 

0.191 

13 

J75 

48469, 13209 

561 

0.45 

1.47 

-0.04 +/- 0.42 

-0.14+/-0.40 

0.410 

0.28 

0.146 

14 

644 

34164, 6255 

954 

1.20 

1.00 






15 

A50 

33386, 9755 

72 

0.60 

2.06 

-0.03 +/- 0.51 

0.05 +/- 0.51 

0.510 

0.24 

0.058 

16 

926 

28578, 8460 

171 

0.81,0.41 

0.87 

0.15+/-0.38 

0.27 +/- 0.39 

0.385 

0.44 

0.309 

17 

461 

28038, 6281 

782 

0.60, 1.02 

1.10 

-0.03 +/- 0.27 

0.14+/-0.27 

0.270 

0.24 

0.143 

18 

A14 

24354, 6448 

115 

0.50 


0.08 +/- 0.41 

-0.06 +/- 0.36 

0.385 


0.100 

19 

J04 

23322, 6460 

188 

1.00 

0.62^ 

0.16+/-0.29 

0.24 +/- 0.30 

0.295 

0.47 

0.288 

20 

A77 

21677, 5423 

318 

0.50 


0.027 +/- 0.63 

0.22 +/- 0.50 

0.565 


0.348 


^Mahabal et al. (2010 . ^Li et al. (19991 ^pry et al. (20121 
(fn) Highlighted columns are CoLiTec users. 


Table 2. The accuracy parameter of thirty observatories that were the most productive in the number of asteroids measurements in 2012 according to the MPC 
data. 


(1) (2) 

N Observatory code 

(3) 

Measurements,objects 

(4) 

Discoveries 

(5) 

D, m 

(6) 

Spix 

(7) 

R.A.(Aa/(TQ,) 

(8) 

Decl.(A5/(T,5) 

(9) 

a" 

(10) 

^pix 

(11) 

ARM 

I 

G96 

2080033, 384204 

17676 

1.50 

1.00 

0.20 +/- 0.33 

0.20 +/- 0.28 

0.305 

0.305 

0.028 

2 

F51 

1948353, 467091 

13785 

1.80 

0.30 

0.07+/-0.15 

0.04+/-0.17 

0.16 

0.53 

0.081 

3 

703 

1723293, 282864 

2278 

0.68 

2.60 

-0.22 +/- 0.65 

0.07 +/- 0.62 

0.635 

0.24 

0.231 

4 

704 

1681504, 262209 

224 

1.00 

2.20 

0.26 +/- 0.67 

0.43 +/- 0.64 

0.655 

0.29 

0.502 

5 

691 

896972, 163714 

7600 

0.90 

1.00 

-0.16+/-0.32 

0.10+/-0.29 

0.305 

0.27 

0.189 

6 

E12 

259295, 62621 

430 

0.50 

1.801 

-0.01 +/-0.51 

0.29 +/- 0.50 

0.505 

0.28 

0.290 

7 

J43 

102641, 22682 

531 

0.50^ 

1.202 

0.19+/-0.48 

0.05 +/- 0.40 

0.44 

0.36 

0.196 

8 

926 

100161, 29986 

454 

0.81,0.41 

0.87 

0.02 +/- 0.37 

0.05 +/- 0.35 

0.36 

0.41 

0.54 

9 

H15 

97878, 24170 

338 

0.45 

2.00 

-0.06 +/- 0.50 

-0.01 +/- 0.53 

0.515 

0.25 

0.061 

10 

106 

72192,17451 

120 

0.60 

2.00 

0.04 +/- 0.36 

-0.12+/-0.34 

0.35 

0.17 

0.126 

11 

A14 

57243, 16239 

159 

0.50 


0.06 +/- 0.37 

-0.02 +/- 0.32 

0.345 


0.063 

12 

J04 

43209,10708 

513 

1.00 

0.62® 

0.21 +/- 0.28 

0.20 +/- 0.27 

0.275 

0.44 

0.29 

13 

DOO 

31494, 7403 

61 

0.40 

2.06 

0,00 +/- 0,57 

-0,06 +/- 0,41 

0,49 

0,23 

0.06 

14 

291 

24272, 6224 

28 

1.80 

0.60 

0.07 +/- 0.33 

0.13+/-0.28 

0.305 

0.50 

0.148 

15 

461 

23847, 5615 

170 

0.60, 1.02 

1.10 

0.00 +/- 0.27 

0.15+/-0.27 

0.27 

0.24 

0.15 

16 

644 

22714,4486 

332 

1.20-* 

1.004 






17 

H21 

22672, 3870 

181 

0.61,0.81,0.76 0.802 

0.03 +/- 0.34 

0.01 +/- 0.36 

0.35 

0.43 

0.032 

18 

141 

21245, 2392 

1790 

1.20® 

1.01® 

0.11 +/-0.23 

-0.03 +/- 0.23 

0.23 

0.22 

0.114 

19 

A24 

18940, 2412 

0 

0.36 

1.40 

0.14+/-0.37 

0.24 +/- 0.33 

0.35 

0.25 

0.278 

22 

A50 

11559,3725 

13 

0.60 

2.07 

0.25 +/- 0.50 

-0.04 +/- 0.46 

0.48 

0.23 

0.253 

iMahabal et al. (20111. 

^Oiv et al. (20121. ^Abreu et al. (20111. ^NEAT-PALOMAR (20131. '^Waszczak et al. (201 31 




(fn) Highlighted columns are CoLiTec users. 









object position in the pixels on the CCD frame (Figures 3b; Col- their class of telescopes (small aperture). In 2011 (2012), this pa- 

umn 10 in Tables 1-3) is used to characterise namely an efficiency rameter was equal to 0.25 (0.25) pixel and 0.24 (0.23) pixel for 

of the mathematical method, which is applied for coordinate 

mea- observatories H15 and A50, respectively, as well as was equal to 

surements. In other words, it allows the observer to be disengaged 0.23 pixel for DOO observatory in 2012. In 2013, the accuracy of 


from parameters of the used CCD matrix and other devices. Ac¬ 
cording to this parameter, the observatory-partners of the CoLiTec 
program have one of the best results among the observatories in 


measurements of the observatories (CoLiTec partners) fell approx¬ 
imately 20% for all the mentioned parameters due to an error in 
software that was immediately fixed and corrected completely. As 
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Table 3. The accuracy parameter of thiity observatories that were the most productive in the number of asteroids measurements in 2013 according to the MPC 
data. 


(1) 

N 

(2) 

Observatory code 

(3) 

Measurements, objects 

(4) 

Discoveries 

(5) 

D, m 

(6) 

Spix 

(7) 

R.A.(Aa/a-Q,) 

(8) 

Decl.(Ai/(Ti) 

(9) 

u" 

(10) 

^pix 

(11) 

ARM 

1 

F51 

2279609, 506894 

14168 

1.80 

0.30 

0.13-t/-0.07 

0.06+/-0.14 

0.135 

0.45 

0.092 

2 

G96 

1950642, 343808 

11908 

1.50 

1.00 

0.04 +/- 0.32 

0.05 +/- 0.28 

0.300 

0.30 

0.064 

3 

703 

1844330, 289086 

1494 

0.68 

2.60 

-0.14+/-0.66 

0.22 +/- 0.64 

0.650 

0.25 

0.260 

4 

691 

742001, 139225 

5594 

0.90 

1.10 

-0.16+/-0.31 

0.12+/-0.30 

0.315 

0.28 

0.200 

5 

D29 

551094, 136964 

262 

1.04 

1.70 

0.03 +/- 0.53 

-0.04 +/- 0.49 

0.510 

0.30 

0.050 

6 

141 

440712, 52579 

2270 

1.201 

l.OH 

0.06+/-0.18 

0.02+/-0.17 

0.175 

0.17 

0.063 

7 

E12 

229747, 48026 

204 

0.50 

1.80^ 

-0.02 +/- 0.50 

0.28 +/- 0.46 

0.480 

0.26 

0.280 

8 

926 

179570, 53662 

750 

0.81, 0.41 

0.87 

0.15+/-0.39 

0.10+/-0.36 

0.375 

0.43 

0.180 

9 

J43 

151983, 27006 

1006 

0.50® 

1.203 

0.11 +/-0.31 

-0.03 +/- 0.29 

0.300 

0.25 

0.114 

10 

W84 

110213,8518 

4160 

403.4 

0.274 

0.13+/-0.13 

0.14+/-0.13 

0.130 

0.48 

0.191 

11 

H15 

107989,25282 

156 

0.40 

2.00 

0.09 +/- 0.62 

0.02 +/- 0.60 

0.610 

0.305 

0.092 

12 

704 

81054, 17833 

4 

1.00 

2.20 

0.29 +/- 0.64 

0.38 +/- 0.63 

0.635 

0.28 

0.478 

13 

J04 

58307, 14670 

576 

1.00 

0.625 

0.25 +/- 0.30 

0.23 +/- 0.28 

0.290 

0.46 

0.340 

14 

DOO 

44658,10850 

34 

0.40 

2.06 

0.01 +/- 0.72 

-0.12 +/- 0.54 

0.630 

0.305 

0.120 

15 

G32 

36416, 4654 

654 

0.40 

1.13 

0.03 +/- 0.35 

0.05 +/- 0.32 

0.335 

0.29 

0.058 

16 

106 

18601, 4502 

67 

0.60 

2.00 

0.04 +/- 0.39 

-0.05 +/- 0.37 

0.370 

0.19 

0.064 

17 

H21 

16924, 2994 

60 

0.61,0.81,0.76 

0.80® 

0.04 +/- 0.33 

-0.04+/-0.31 

0.320 

0.4 

0.002 

18 

461 

15688, 3787 

110 

0.60, 1.02 

1.10 

-0.02 +/- 0.24 

0.17+/-0.27 

0.255 

0.23 

0.171 

19 

644 

15221, 3317 

63 

1.20i 

i.ooi 






20 

291 

15197,4002 

1 

1.80 

0.60 

0.02 +/- 0.35 

0.14+/-0.32 

0.335 

0.55 

0.141 


iWaszczak e.t al. (•2013) . ^Mahabal et al. 12011) . ^Qry et al. (2012) . ^Honscheid et al. (2008) . ^Abreu et al. (2011) . ^Li et al. (1999) . ^NEAT-PALOMAR (2013) 
(fn) Higbligbted columns are CoLiTec users. 


a result, the observatory-partners of the CoLiTec program returned 
to their previous positions as for the indexes of measurement accu¬ 
racy (see, in detail, the current MPC site). It is important to empha¬ 
sise that standard deviation api^ of object position (Figures 3b) is 
an artificial parameter. It is not a strong objective because such fac¬ 
tors as the exposure time, the telescope optical scheme, the height 
above sea level, and many others are not taken into account in it. For 
example, according to this parameter, the Pan-STARRS 1 observa¬ 
tory (F51) had lost positions, although it had the best astrometry 
accuracy among all the asteroid surveys. We may suggest that an 
implementation of our method, which is free from a possible loss 
of measurement information contained on the CCD frames, when 
using at this observatory most likely can yield the best results than 
that was used. 

We also compared the CoLiTec software and the Astrometrica 
that is widely accepted among amateurs. 

Both the softwares: realise such functions as frame calibra¬ 
tion; support of astrometric and photometric catalogues of stars 
in the local and online modes; record the coordinate information 
(WCS) in the FITS-frame title; have interactive mode for the ob¬ 
ject position measurements; have a magnifier tool; have an auto¬ 
mated search of moving objects; have a mode for visual inspection 
of detected moving objects; provide output of astrometric measure¬ 
ments in the MPC format; transfer data mode from the interface to 
the MPC; display the known and discovered objects on the CCD 
frame. However, the CoLiTec does not realise such functions as 
Track&Stack for adding the frame, and identification of detected 
moving objects with a local database MPCORB. 

At the same time, unlike the Astrometrica, the CoLiTec: re¬ 
alises such functions as astrometric reduction of the CCD frames 
with large fields of view (2‘^®® or more); separates application of 
the reference catalogues for astrometric and photometric reduc¬ 
tions; provides an automated search of moving faint objects (SNR 
2.5); has automatic data processing; saves results on the process¬ 


ing frames (the real detected objects; objects, which are rejected by 
the operator, etc.); identifies the detected moving objects with the 
online database MPCORB (MPC); identifies the stationary objects 
with database of variable stars (VSX) and galaxies (HyperLeda); 
has a software modular design (the ability to connect the individ¬ 
ual modules). The CoLiTec software gives more accurate measure¬ 
ments of faint celestial objects as well as contain a more reliable 
method for identifying the frames with a reference star catalogue 
allowing to improve (sometimes significantly) an accuracy of the 
object position. 

Along with the analysis of indexes of measurement accuracy 
(see, the current MPC site), we conducted a comparative analysis of 
the accuracy of both the softwares (CoLiTec vs. Astrometrica) af¬ 
ter processing the same frame. We selected 19 series for each of the 
4 frames. A preliminary analysis included 36 series. However, the 
rest of the series had no reliable identification with the used star cat¬ 
alogue. Also, we excluded frames with significant disruptions in the 
daily maintenance and frames taken at very high wind. All frames 
were obtained at the observatory ISON-NM (H15) with help of a 
40-cm telescope Santel-400AN and CCD-matrix FLIML09000-65 
(3056 X 3056 pixels, a pixel size of 12 microns) in the period from 
March 4, 2014 to March 30, 2014. The exposure time was 150 sec¬ 
onds. 

We used only such object positions, the m easurements o f 
which are included in the MPCAT-OBS archive dMPCAT-OBSb . 
The measurements, however, were reprocessed with the CoLiTec 
software. To set the reference values of object positions on the 
time of measurem ents formation, we used the HORIZONS ser¬ 
vice i IHORIZONSI) . In total, we used 2002 measurements (mea¬ 
surements with CoLiTec were done on 253 more). The Astromet¬ 
rica in several cases issued a sign about impossibility to guarantee a 
reliable measurement of object position (Centroid = -1). It is often 
associated with an attempt to measure the position of star trails or 
faint objects involved with a brighter star. More than half of them 
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Figure 3. The comparative analysis of the accuracy parameters of object position for the most productive observatories in the number of asteroids measurements 
in 2011 (left), 2012 (middle), and 2013 (right) by: a) module of the average residuals of object position measurements; b) standard deviation estimations o-pi^ 
of object position, in pixels; c) standard deviation estimations of object position in the right ascension, cia, in arcseconds; d) standard deviation estimations of 
object position in the declination, as, in arcseconds. 


have an SNR not exceeding 3.5. The results of a comparative anal¬ 
ysis are shown in Table 4. One can see that measurements with the 
Astrometrica at a low SNR have an RMS of 30-50% larger than 
that of the CoLiTec (see, also Fig. 4 and 5). Mean deviations in 
measurements with the CoLiTec and the Astrometrica are the same 
in general, and so these data are not sh own. The more detailed co m- 
parative analysis is given in our paper dSavanevveh et al.l l l2015h ). 


5 CONCLUSIONS 

A new iteration method of asteroid coordinate estimation on the 
digital image has been developed. The method operates by con¬ 
tinuous parameters (asteroid coordinates) in a discrete observation 
space (the set of pixels potential of the CCD matrix). 

High rates of the CoLiTec program during 2011-2012 as con¬ 
cerns with the accuracy of measurements have been obtained due 
to the use of the subpixel Gaussian model. This model of the ob¬ 
ject image takes into account a prior form of the object image and 


Table 4. Comparative analysis of the softwares CoLiTec and Astrometrica 
as concerns with measurements of object positions for the numbered aster¬ 
oids 


(1) 

Deviation 

(2) 

Astrometrica 

(3) 

CoLiTec 

The average deviation of RA ,'' 

0,11 

0,11 

The average deviation of DE, " 

-0,04 

-0,03 

RMS deviation of RA, " 

0,77 

0,50 

RMS deviation of RA, " 

0,67 

0,39 


consequently it is adapted more easily to any forms of real image. 
In other words, nevertheless that a real coordinate distribution of 
photons hitting the pixels on the CCD frame is not known, the 
form of this distribution is known a priori and its parameters can 
be estimated according to the real object image. At that time many 
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Figure 4. Distribution of deviations of the equatorial coordinates of object 
by SNR ranges (Astrometrica vs. CoLiTec) 
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Figure 5. Distribution of deviations of the equatorial coordinates of object 
by magnitude ranges (Astrometrica vs. CoLiTec) 


other methods mentioned in Introduction consider by default that 
the density of hit photons inside the pixel is uniform. 

The advantages of subpixel Gaussian model become more 
obvious for the fainter celestial objects. Moreover, the developed 
method has a high measurement accuracy as ■well as a lo'w cal¬ 
culating complexity because the maximum likelihood procedure 
is implemented to obtain the best fit instead of the least-squares 
method and Levenberg-Marquardt algorithm for the minimisation 
of the quadratic form. 

The efficiency of the proposed method, including its advan¬ 
tages for accurately estimating asteroids coordinates, was con¬ 


firmed during observations as the part of the CoLiTec (Collection 
Light Technology) program for automatic discoveries of asteroids 
and comets on a set of the CCD frames. Efficiency is a crucial fac¬ 
tor in the discovery of near-Earth asteroids (NEA) and potentially- 
hazardous asteroids. Current asteroid surveys yield many images 
per night. It is no longer possible for Ihe observer to quickly view 
these images in the blinking mode. It causes a serious difficulty for 
large-aperture wide-held telescopes, capturing up to several tens of 
asteroids in one image. The CoLiTec software solves the problem 
of the frame processing for asteroid surveys in a real-time mode. 
We compared also our software and Astrometrica, which is widely 
used for detecting the newly objects. Limits of measurements of 
CoLiTec software are wider than of Astrometrica, but that the most 
valuable, this expansion comes into area of extremely small SNR 
allowing to search the fainter Solar System small bodies (measure¬ 
ments with the Astrometrica at a low SNR have an RMS of 30- 
50% larger than that of the CoLiTec). As for the area of SNR>7, 
the results of CoLiTec and Astrometrica are approximately iden¬ 
tical. However, exactly the area of extremely small SNR is more 
promising for the discovery of new celestial objects. 

The automatically detected small Solar System bodies are sub¬ 
ject to follow-up visual conhrmation. The CoLiTec software is in 
use for the automated detection of asteroids at the Andrushivka As¬ 
tronomical Observatory, Ukraine (since 2010), at the Russian re¬ 
mote observatory ISON-NM (Mayhill, New Mexico, USA) since 
2010, at the observatory ISON-Kislovodsk since 2012, and at the 
ISON-Ussuriysk observatory since 2013 (se e Tables 1-3 ). As th e 
result, four comets (C/2010 XI (Elenin) l Elenin_et_aL 
P/2011 NOKElenin) (lEleninetal.l ( 1201 ih 


Elenin et al. 


I201C1). 

1 20131) ). 


C/2012 SIHSON) jNevskil (|2012|) ). and P/2013 V3 (Nevski) 

In evskil (I 2 OI 3 I) ) as well as more than 1500 small Solar System bod¬ 
ies (including live NEOs, 21 Trojans of Jupiter, and one Centaur 
object) have been discovered. 

In 2014 the CoLiTec software was recommended to all mem¬ 
bers of the Gaia-EUN-SSO network (a network for Solar Sys¬ 
tem transient Objects) for analysing observations as a tool to de¬ 
tect faint moving objects in frames. Information about CoLiTec 
with link to web-site has been posted on the Gaia-EUN-SSO Wiki 
(https : //www. imcce . fr/gaia-fun-sso/l. 

The authors thank Dr. E. Velichko for his useful comments. 

We are grateful to the reviewer for his helpful remarks that 
improved this work. We express our gratitude to Mr. W. Thuillot, 
coordinator of the Gaia-EUN-SSO network, for the approval 
of CoLiTec as a well-adapted software to the Gaia-EUN-SSO 
conditions of observation. We also thank Dr. Ya. Yatskiv for 
his support of this work in frames of the Target Program of 
Space Science Research of the National Academy of Science 
of Ukraine (2012-2016) and the Ukrainian Virtual Observatory 
(http://www.ukr-vo.orgl. The CoLiTec program is avail¬ 
able through http://colitec.neoastrosoft.com/en 
(one can access to the download package 
http://www.neoastrosoft.com/download_en and 
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